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Abstract 

The paper presents an analysis of the time reversal in independent-multipath Rayleigh- 
fading channels with N inputs (transmitters) and M outputs (receivers). The main 
issues addressed are the condition of statistical stability, the rate of information 
transfer and the effect of pinholes. The stability condition is proved to be MC <C 
N c f[B for broadband channels and M <C N c g for narrowband channels where C 
is the symbol rate, B is the bandwidth and N c r is the effective number (maybe 
less than 1) of transmitters. It is shown that when the number of screens, n — 1, 
is relatively low compared to the logarithm of numbers of pinholes N e s is given 
by the harmonic (or inverse) sum of the number of transmitters and the numbers 
of pinholes at all screens. The novel idea of the effective number of time reversal 
array (TRA) elements is introduced to derive the stability condition and estimate 
the channel capacity in the presence of multi-screen pinholes. The information rate, 
under the constraints of the noise power v per unit frequency and the average total 
power P, attains the supremum P/v in the regime M A N c r 3> Pj {yTS). In particu- 
lar, when iV e fj ^> M S> P/{Bv) the optimal information rate can be achieved with 
statistically stable, sharply focused signals. 
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Fig. 1. MIMO Broadcast Channel 

1 Introduction 



Time reversal (TR) of waves has received great attention in recent years and 
been extensively studied for electromagnetic [2] , [20] as well as acoustic prop- 
agation (see [14] and the references therein). A striking effect of time reversal 
in randomly inhomogeneous media is the superresolution of refocal signals 
[1], [13] which implies low probability of intercept and holds high potential in 
technological applications such as communications [9], [25], [10], [17], [18]. 

An issue prior to superresolution, however, is statistical stability, namely the 
question: How many antennas and how much bandwidth does one need to 
achieve self-averaging in TR so that the received signals are nearly deter- 
ministic, independent of the channel statistics? In this paper we answer this 
question for independent-multipath Rayleigh fading channels, with multiple 
inputs and multiple outputs (MIMO), commonly used in wireless communi- 
cation literature, see, e.g. [23]. We also introduce the novel idea of effective 
number of transmitters to analyze the effect of multi-screen pinholes on sta- 
bility and capacity. 

In the MIMO-TR communication scheme as studied in [9], [12], the M well- 
separated receivers first send a pilot signal to the iV-element time reversal 
array (TRA) which then uses the time-reversed version of the received signals 
to modulate the data symbols and retransmit them back to the receivers. 
One of the main results obtained here is that the time reversal process is 
statistically stable when 
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MC <C N e gB, for broadband channels 
M <C N cS , for narrowband channels 



(1) 
(2) 



where C(< 2B) is the symbol rate, B is the bandwidth and N eS is the effective 
number of transmitters (maybe less than one). In the presence of (n — l)-screen 
pinholes, we show that the effective number of transmitters is asymptotically 
the harmonic sum of the number of transmitters and the number of pinholes 
of every screen when all these numbers are greater than 2™. 

The LHS of (1) is the number of degrees of freedom per unit time in the 
constellation of input data-streams while the RHS of (1) is roughly the number 
of degrees of freedom per unit time in the channel state information (CSI) 
received by TRA from the pilot signals. The latter has to be larger than 
the former in order to reverse the random scrambling by the channel and 
achieve deterministic outputs. The stability condition iV 3> 1 for narrow- 
band channels or B ^> f} c (the coherence bandwidth) for broadband channels, 
when M is small and the pinholes are absent, have been previously discussed 



In Section 4 and 5.2, we analyze the information rate of the TR broadcast 
channel in the presence of noise. We show that the optimal information rate 
R ~ P/v, under the power and noise constraints, can be achieved in the regime 
M A TVeff 3> Pj iyB) where v is the magnitude of noise per unit frequency and 
P the average total power input. In particular, when iV cff 3> M 3> Pj{Bv) 
the optimal information rate can be achieved with statistically stable, sharply 
focused signals. 



2 TR-MIMO communication 

First let us review the MIMO-TR communication scheme as described in [12] 
which is an example of broadcast channel [23]. 

The M receivers located at y 3 -, j = 1, M first send a pilot signal / e tult g{B^ 1 (uo— 
uj ))duj5(x — yj) to the iV-element TRA located at Xj,i = 1,...,N which 
then uses the time-reversed version of the received signal / e luJt g(B~ l (u — 
uo))H(yj,Xi;uj)duj to encode a stream of symbols and retransmit them back 
to the receivers. Here H is the transfer function of the propagation channel at 
the frequency uo from point y to x and g 2 (uo) is the power density at uo. Let 
H(a;) = [Hij(uj)], Hij(uj) = -f/"(xj, y,-; u), be the transfer matrix between the 
transmitters and receivers. The reciprocity implies that H(yj,Xi;u) = Hij(u) 
and H*(uj) = H(— uo) where * stands for complex conjugation. Let us assume 
that g is a smooth and rapidly decaying function such as the Gaussian func- 
tion. Naturally the relative bandwidth B/u is less than unity so that u ^> 1 



in [1], [7], [8], [9], [20]. 



3 



if B 1. In this paper we will assume B/ujq <C 1 to simplify the frequency 
coherence structure below (Section 3). We have chosen the time unit such that 
the speed of propagation is one and the wavenumber equals the frequency. 

The signal vector S = (Sj) arriving at the receivers with delay L + 1 is then 
given by [12] (see also [1], [6]) 



W M „ N 



^W = EE m »W / e-^^(^^)E^H^H^ (3) 



1=1 i=l " k=l 



where mj(ri),l = 1,...,W < oo is a stream of symbols intended for the j- 
th receiver transmitted at times T\ = It, t > 0. In vector notation, we have 
S = Ei=i/e-^*- Ti ^( J B- 1 (cj-cj ))HHt(cu)m(r i )^ where is the conjugate 
transpose of H and 111(77) = ( m j( T i))- Let us note that while all the TRA- 
elements are coordinated and synchronized the receivers do not know the chan- 
nel and can not coordinate in decoding the total signal vector received. As a 
consequence, the multi-user interference arises and can be a serious impedance 
to communications. An advantage of the time reversal scheme is the possibility 
to use the (statistical) stability property to achieve the following asymptotic 



_ N 

- Mt - n) 9(^^) E H^H^duj ~ BSve-^^r-^KBh - t)) 

k=l 



so that Sj(t) ~ BYZi^^^'^'^^laKBin ~ *)) and each receiver 
receives the input symbols with little interference. Here and below JF _1 stands 
for the inverse Fourier transform. 



3 Statistical stability 



One of the main goals of the present note is to characterize the stability regime 
for the independent-multipath Rayleigh fading channel in which Hij(u) are in- 
dependent CAf(Q, <r), the zero-mean, variance-a circularly symmetric complex- 
Gaussian random variables and {Hij{uo)} i ■ u are a jointly Gaussian process. 
The independent-multipath Rayleigh fading is an idealized model for richly 
scattering environment, after proper normalization, when the spacings within 
the transmitters and receivers are larger than the coherence length t c of the 
channel. In general, the coherence length is inversely proportional to the angu- 
lar spread [23] and sometimes can be computed explicitly in terms of physical 
properties of the channel [12] . For diffuse waves the coherence length is known 
to be on the scale of wavelength [28], [27]. 
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We set the variance a = l/( N V M) so that the average input power is no 
less than the average output power. The value of a would not change the 
conditions of statistical stability but will affect the discussion of information 
transfer in the next section. 

Let us calculate the mean and the variance of the signals with respect to 
the ensemble of the channel. Let E denote the channel ensemble average. For 
simplicity, we assume that \mi{ji)\ = By the Gaussian rule for the 

calculation of moments we have 

w 

ES = BNamY,^ o{t - Tl) J r - l [g}{B{t-T l )). (4) 
1=1 

Let r > (2B)- 1 so that the summation in ES is _B-uniformly bounded as 
W -> oo. 

The statistical stability of the signals can be measured by the normalized 
variance of the signals at the receiver j 

V,-(r B ) = |E |^ r V 3 {r n ) = E\S s \ 2 (r n ) - |E^-(r n )| 2 , 

Vj, n and the time-reversed signals are stable when V,(r n ) — > 0, Vj, n. Note 
that VJ 1 (ti) is exactly the signal-to-interference ratio (SIR) at receiver j. 

Let f3 c be the coherence bandwidth of the channel such that 
E[/^(u;)#^,(u/)] w <7/(w , ^-^-)5«'^y 

where f(u>o, •) is a continuous, rapidly decaying function and /(c^o, 0) = 1 (see 
[11], [12] for a rigorous example). Here we have used the fact that the relative 
bandwidth B/cu is small so that / is independent of the precise value of 
the frequency. Below we shall suppress the argument cu in /. The coherence 
bandwidth f3 c is inversely proportional to the delay spread and hence the 
delay-spread-bandwidth product (DSB) is roughly Bp- 1 [11], [12], [23]. In the 
diffusion approximation f3 c is given by the Thouless frequency DbL~ 2 where 
Db is the Boltzmann diffusion constant, equal to the energy transport velocity 
times the transport mean free path, and L the distance of propagation [19], 
[29]. 

The broadband, frequency-selective (BBFS) channel is naturally defined as hav- 
ing a large DSB, i.e. Bfc 1 > I. Since B < uj , uj e [uj - B/2,uj + 5/2] and 
— uj are separated by more than [3 C in a BBFS channel. On the other hand, 
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B <C p c corresponds to the narrow-band, frequency-flat (NBFF) channel. For 
convenience in the subsequent analysis, we shall think of the NBFF channel 
as the limit f3 c — > oo and the BBFS channel as the limit f3 c — > while u , B 
are fixed. In either case, we have 



VG-(0»^ 2 E E m^mtin.y^-^ (5) 

i=l Z,Z'=1 

Consider the NBFF channels first. We obtain by passing to the limit j3 c — > oo 



tyf) « iW 2 B 2 |/| 2 (0) E | E^K^^W^ - tj)) 

i=i z=i 

In view of (4) the stability condition N ^> M for NBFF channels then follows 
easily. On the other hand, the BBFS channels (f3 c — > 0) yields 



m w 

i=l /,/'=! 



M W 

« 7V(7 2 fl& E E ™^)^ -1 [l/| 2 ] (A(rZ - 0) 
i=i ;=i 

w 

x E;<(^)c iwor(, - , ' ) ^ 1 [bl 2 ](Sr(Z - I')). (6) 
Z'=l 

Several observations are in order. First, due to r > (2B) _1 the summation over 
/' in (6) is convergent as W — > oo uniformly in £?. Second, due to the term 
.F _1 [|/| 2 ](/3 c (tZ — i)), there are effectively C/3" 1 terms in the summation over 
/ where C = r~ l is the number of symbols per unit time in each data-stream. 
As a result, the variance Vj^N a 2 BMC fl 2 is independent of f3 c . It then follows 
that Vj — > if and only if NB ^> MC for BBFS channels. The transition 
to the condition N ^> M for NBFF channels takes place when B ~ C, i.e. 
t~B-\ 



The stability condition can be interpreted as follows: NB is the number of 
degrees of freedom in the CSI collected at the TRA per unit time; MC is 
the number of degrees of freedom in the ensemble of messages per unit time; 
the stability condition NB 3> MC says that in order to recover the input 
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messages, independent of the channel realization, and thus reverse the random 
scrambling by the channel, the former must be much larger than the latter. 
In light of this interpretation, the stability condition derived above appears to 
be sharp. 

A detailed, rigorous analysis of the MIMO-TR channel modeled by a stochastic 
Schrodinger equation, in the parabolic approximation of scalar waves, with a 
random potential is given in [12]. 



4 Rate of information transfer 

In this section we discuss the information rate for a memoryless channel which 
is constructed out of the time-invariant channel model analyzed in Section 3. 
The temporal dependence is introduced by drawing an independent realization 
from the Rayleigh-fading ensemble of transfer matrices after each use of the 
channel, i.e. after each delay spread (or two if the time for channel estimation is 
included) . This is obviously an idealization but widely used in communications 
literature [31], [15]. The coherence time of the resulting ergodic channel model 
is then much longer than one delay spread. We assume as in standard practice 
that in addition to the random channel fluctuations additive-white-Gaussian- 
noise (AWGN) is present at each receiver, that the input signal vector is 
multivariate Gaussian and that the channel, the noise and the input signal 
are mutually independent. 

For the Rayleigh fading channel prior to adding noise, each frequency compo- 
nent of the time reversed signal Sj in (3) 

M N 

i=i fc=i n 

N _ N _ 

1 - itU^l Z • 

iV-degree central x 2 r - v - N(M - 1) i.i.d. zero-mean r.v.s 

is a sum of a central \ 2 random variable with iV degrees of freedom and 
N(M—1) i.i.d. mean-zero random variables. This is due to the assumption that 
different entries of the transfer matrix are mutually independent zero-mean 
Gaussian random variables. Therefore, for N ^> 1 the interference statistic is 
approximately Gaussian, by the Central Limit Theorem. More generally, after 
synthesizing all the available frequencies, the interference statistic becomes 
approximately Gaussian if NBfl^ 1 ^> 1 which is always the case for the BBFS 
channels. In a BBFS (resp. NBFF) channel, NB(3~ l (resp. N) is the number 
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of independent subchannels from TRA to each receiver. 



Moreover, each frequency component of Sj has the mean 



M N 

E [EE 



, UJ — OJq , 



/, I B 



UJ — Uq s 



mi(n)g( — ^-)H jk (uj)H* k (uj) = Nag(— ^)m,(r,). (7) 



which exhibits the simple input-output relation: The w-component of the in- 
put signal for the j-th receiver is rrijg(uj) and the received signal component 
is Namjg(uj) corrupted by the noise and interference which for N ^> 1 is 
approximately Gaussian. Since the M receivers operate independently of one 
another, the total time-reversal broadcast channel consists of M independent 
subchannels in parallel each of which has the above input-output relation. 
Thus the total information rate is the sum of those of the M subchannels 
from TRA to individual receivers. And, in view of the simple input-output re- 
lation, each subchannel can be viewed as a single-input-single-output (SISO) 
linear filter channel corrupted by (approximately) Gaussian noise/interference 
for which Shannon's theorem is applicable. 

According to Shannon's theorem [5] the ergodic capacity (in nats per unit time 
and frequency) of a SISO linear filter channel is In (1 + SINR) where SINR, the 
signal-to-interference-and-noise ratio at each receiver, is given by the harmonic 
sum of the SIR, the signal-to- interference ratio and SNR, the signal-to-noise 
ratio, i.e. SINR = (SIR -1 + SNR -1 ) -1 . For extension of Shannon's result to the 
MIMO setting, see [15], [31]. 

Analogous to the NBFF channels in Section 3, SIR(w) = Vj 1 ~ N/M, inde- 
pendent of /i and uj. Let v be the noise power, per unit frequency, at each 
receiver. Suppose the average transmission power is constrained to P and all 
the transmit and receive antennas are identical. 

Since the value of a would affect SNR (but not SIR) we discuss the two cases 
N > M and N < M separately. 

Case 1: N > M. In this case, a = N" 1 and in view of (7), SNR(cj) = /i 2 /V 
where \x = \mj\ can be related to the total power constraint P as fi 2 M ~ P/B 
since the average input power per unit frequency is 

N M 

E E WtOI 2 !*?! 2 ^- 1 ^ - uj ))E\H lk (cu)\ 2 ~ MNafi 2 = M/i 2 . 

k=l i=l 



Thus SNR(c<j) ~ P/(uBM). Therefore the total channel capacity (in nats per 
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unit time) is roughly given by 



(8) 



Now we ask the question: What is the maximal rate at which a TRA, with 
fixed number of elements N, fixed average total power P and fixed noise level 
(per frequency) v, can transfer information if there is no limitation to the 
number of receivers M and the bandwidth Bl 

Expression (8) can be optimized at the limit M»P/ {vB) to yield the optimal 
information rate of P/v which is linearly proportional to the power. We see 
that the simplest strategy for optimizing the information rate of a given TRA 
under the the power and noise constraints is to enlarge the bandwidth B as 
much as possible. And if we can satisfy N 3> M 3> P/(vB) then we can 
achieve stability as well as the optimal information rate. 

Consider the thermal noise power v = k B T where k B is the Boltzmann's 
constant and T the temperature. Then the above result implies that the energy 
cost per nat is P/R ~ k B T which is consistent with the classical result of 
minimum energy k B T requirement for transmitting one nat information at 
temperature T [24], [21]. 

Case 2. N < M. In this case, a = M' 1 and (7) implies that SNR ~ 
N 2 fx 2 /(M 2 v) where n is related to P by fi 2 = P/{NB). Hence SNR ~ NP/{M 2 vB). 
With SIR ~ N/M and Shannon's theorem, the channel capacity is roughly 



which achieves the optimal rate Pj v in the regime N = M P/ (Bis). On the 
other hand, for M <C Pj (Bv), the information rate becomes BMln (1 + N/M) < 
BN which is much smaller than Pj v. 

Therefore we conclude that under the power and noise constraints the condi- 
tion for the optimal information rate P/v is iV > M ^> P/(Bv), which can 
be achieved by sufficiently large bandwidth, whereas the additional condition 
iV ^> M, which, sufficient for the Gaussian approximation to the interference 
statistic, would also guarantee stability. 

Before ending this section, let us compare the capacity in the conventional, 
non-TR MIMO channel as calculated in [15], [31], [22], [30]. Consider the 
non-TR single-user channel with the M transmit antennas (on the right of 
Fig. 1) which have no channel knowledge and the iV(> M) receive antennas 
(on the left of Fig. 1) as the single user which has perfect CSI. This is, of 




(9) 
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Fig. 2. Single-screen pinholes 

course, the reciprocal case of the TR broadcast channel. In this case, SNR ~ 
Pj{MBv) and it is shown in [15] and [31] that the ergodic capacity of the 
single-user narrowband Rayleigh-fading channel scales like .BMlnSNR at high 
SNR which can be recovered from (8) by imposing the additional constraint 
M < P/ {vB) < N. And as we learn from the discussion of Case 1 above, this 
is not the regime for achieving the optimal information rate Pj v. 

The same results as discussed in this section are obtained for the parabolic 
Markovian channel model in [12]. 



5 Pinhole effect 

Pinholes are degenerate channels that can occur in a wide family of channels, 
outdoor as well as indoor, see Fig. 2 and 3. While preserving the co-channel 
decorrelation, pinholes have been shown to severely limit the degrees of free- 
dom and reduce the channel capacity [3], [16], [4]. In this section, we introduce 
the notion of effective number of TRA elements to analyze the multi-screen 
pinhole effect on TR in Rayleigh fading. 

Let us begin with the simplest case of single-screen pinholes as illustrated 
in Fig. 2. Let h^(u) be the N x K transfer matrix from the TRA to the 
pinholes and h( 2 )(u;) the K x M transfer matrix from the pinhole to the M 
receivers at frequency uj. The combined channel can be described by U(u>) = 

hW(w)hW(w) = [£^A-fcM4jM] in which 4?H and /igV) are as- 
sumed to be independent CAf(0, <7i) and CAf(0, a 2 ), respectively, and {h^\uj), h\f(uj)} i ■ 
are jointly Gaussian processes. To prevent the average input power from being 
less than the average output power we set IE | j | 2 = Ko\Oi = (iVV M) _1 , Vi, j . 
Note that the entries of H are in general not independent r.v.s. 

As before we assume the frequency coherence structure 

E[h%\u,)h$>(u;')}*a k f(^^)5 u ,5 jf , Mk (10) 
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where, for simplicity, / and (3 C are taken to be independent of the screens. 
Straightforward calculations with the Gaussian rule show that the mean signal 
is 

w 

E[Sj(t)] = BNK<T 1 a 2 y £m j ('n)r- 1 \9](B('n - t)) 

1=1 



and its variance is 



i,v=i J B B p c 

M —(/ M 

x (m^m^rv) + N^m^m*^,) + K\f\\^-) £ m^m*^)) (11) 



i=i 



In view of the observations following eq. (6) we have the estimate Vj(t) ~ 
B 2 KN(MN+MK+l)ala%\n\ 2 for the NBFF channels and Vj(t) ~ BCKN(MN+ 
MK + l)cr^crf |/x| 2 for the BBFS channels. As in (6) the variance does not de- 
pend on the coherence bandwidth f3 c . Therefore we obtain the normalized 
variance of the signal to the leading order (N, A"))!) 



| M(iV -1 + K~ l ), for the NBFF channels 

| MCB-^N- 1 + K- 1 ), for the BBFS channels. 



The result suggests the notion of effective number of TRA-elements given 
by N e g = NK(N + K)~ l , namely the harmonic sum of N and K, so that 
Vj ~ MCB^Nj for the BBFS channels and Vj w MN~g for the NBFF 
channels. For TV, K ^> 1 the number of statistically independent paths is 
roughly iV eff x M. 

The previous case without pinholes corresponds to the limiting case K ^> N. 
For a fixed K, however, the previous benefit of stability with large number 
of TRA elements (N ^> 1) disappears. The multiple antennas in TRA are 
essentially screened out by the pinholes and the effective number of TRA- 
elements becomes K. 



5. 1 Multi-screen pinholes 



The same analysis can be applied to channels with (n — 1) screens of pinholes 
such as illustrated in Fig. 3. Let Kk, k — 1, ..n — 1 be the number of k-th screen 
pinholes. Let be the transfer matrix for the k-th stage channel whose 
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Fig. 3. Multi-screen pinholes 

entries are independent CA/"(0, <Jk) and let the transfer matrices of different 
stages be mutually independent. Again, in order for the average input power 
to be no less than the average output power we set 

E|/y 2 = K 1 --- K n _ x a x ■ ■ ■ a n = (N V M)~\ (12) 



The condition of statistical stability, however, is independent of the values of 
a k , k = 1, ...,n. 

As noted previously the the normalized variance does not depend on (5 C and its 
order of magnitude is determined solely by the same- frequency moments which 
will be the focus of the subsequent calculation. The calculation of the mean 
is straightforward: E(Htfm)j = NKi • • • K n ^\<j\ • • • <J n mj. Let us analyze the 
second moment of entry a 



EfHHtm) (HH+m)* 

= El V h {n) h {n - l) ■ ■ ■ h m - h {1) - h (1) * h {2) * ■ ■ ■ h {n ~ l) * h (n) * ■ m ■ 

^\ ' l ai„ n i n ,l n - 1 lb l3,l2 lb t2,ll lb nM lb 3Z,32 U 3n,3n-l"'3n+l,3n'"'3 



»l>"'*n 



X 



£ h^h^P* ■ ■ • 4 2 > h$» hf, hf 3 , ■ ■ ■ hf-*> hf . m*, }. 

t < «t n l n' l n -l l 3'2 2' 1 ->2 ' 1 -'3'-'2 JnJ„-l J n + \>Jn J n+1 ) 



According to the Gaussian rule for computing moments, the above expression 
can be represented by 2 n diagrams of 4n vertexes and 2n edges. We distinguish 
two types of edges: the arcs, connecting (un)primed indices to (un)primed 
indices, and the ladders, connecting unprimed indices to primed indices, see 
Fig. 4. 

When a new screen of pinholes, represented by h^ n+1 \ is added, the number 
of diagrams is doubled: one half of them contain the ladders connecting h^l\ 

to h^ 1 ^* and h^^* , , to h^ +1 \, while the other half contain the arcs con- 

at n+l 3n+2,3n+l 3 n +2<3n+l 
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n n-1 / / s n \ \ (n-1)* n* 

• • v V • • • » • • 

n-2 2 1 1* 2* (n-2)* 
(n-2)* 2* 1* 1 2 n-2 

• • • • • a 4 • 

n* (n-1)* \v \ V J y / n-1 n 

Fig. 4. Separable diagram: * means complex conjugation; the top indices are un- 
primed and the bottom indices are primed. 

necting h^ +1 } to h^ + ^* , , and h^ 1 ^* to h^ +1 \, . Straightforward calculation 

with (10) yields the following rule: A new pair of arcs add to diagrams with 
outermost arcs the K 2 (multiplicative) factor and diagrams with outermost 
ladders the K n /M factor; on the other hand, a new pair of ladders add to 
diagrams with outermost ladders the K\ factor and diagrams with outermost 
arcs the K n M factor. 

That is, the diagrams that correspond to the highest power in K±, K 2 , ■ ■ •, 
have the least number of edge-type alternating. Hence for K\, ■ ■ ■ , K n -\ 3> 
2™ ^> N the leading order term in the variance corresponds to the diagram 
with all ladders and is of order K\ ■ ■ ■ K^_ 1 NM while the square of the mean 
corresponds to the diagram with all arcs and is of order K\ ■ ■ ■ K^^N 2 . The 
stability condition thus remains the same as in the case without pinholes. 

Let us consider the more interesting regime in which N, Ki, .., K n _i ^> 2™. We 
claim that to the leading order the normalized variance of the signal is given 
by Vj ~ MCB~ l N~ff where the effective number of TRA-element N e g is given 
by 



namely the harmonic sum of N, Ki, ■ ■ ■ , K n -i. We sketch the proof here. The 
leading order terms in the variance after expectation correspond to the separa- 
ble diagrams in which the arcs are nested and are flanked by the ladders, Fig. 
4. Except for the diagram with all ladders, the separable diagrams all have 
the innermost arcs connecting to h^J^ and hfjP*, to h^\, , which give rise 
to the factors iV 2 (an extra iV than otherwise), and, except for the diagrams 
with all ladders or all arcs, the separable diagrams change the edge-type ex- 
actly once (from arc to ladder). When iV is comparable to K 1: • • • , K n _ 1: the 
contributions from the separable diagrams are comparable to that from the 
diagram of all edges. 

Collecting the terms corresponding to the separable diagrams we have 
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Yl Tl — 1 TL — 1 71 — 1 

/j, 2 NM f[^Y[K j (Y[K k + Nj2K 1 ---K i --- K n - X ) 

i=l j=l k=l i=l 



where Ki means that Ki is absent in the product. Dividing it by iV 2 IIT^i 1 
and accounting for the temporal aspect of transmission as in the observations 
following eq. (6) we obtain the claimed result. 



5.2 Information rate with pinholes 



The notion of the effective number of TRA elements is useful in estimating the 
channel capacity as well as the stability condition in the presence of pinholes 
since SIR is given by N cS /M with C = 2B. 

As the (spatial) subchannel from TRA to each receiver is the sum of NK X K 2 ■ ■ ■ K n _i 
paths which are not necessarily independent, the simplest way for realizing 
Gaussian interference statistic is to assume large degrees of freedom in fre- 
quency BP' 1 1 so that each spatial subchannel gives rise to a sum of Bfi~ x 
roughly i.i.d. r.v.s. This works only for the BBFS channels. For the NBFF 
channels, we assume the worst-case scenario if m i n = min [N, K — 1, • • • , K n _i\ ^> 
1 because each subchannel can be regrouped into a sum of NKiK 2 ■ ■ ■ K n ^i/ K min 
terms each of which is a sum of K min i.i.d. r.v.s. 

Due to the normalization (12) the input-output relation in (7) and the dis- 
cussion in Section 4 (Case 1 & 2) remain valid if iV is replaced by N cS . In 
particular, the same optimal information rate P/u is achieved in the regime 
N cS AM^>P/(Bu). 

As analyzed before, when the condition N, Ki, .., K n _ 1 ^> 2 n is satisfied, N c g 
is the harmonic sum of N,K\, ...,K n _i and therefore we have the estimates: 
Kmin/n < N c g < K max /n where K m { n and K max are the minimum and maxi- 
mum of N, Ki, K n _i, respectively. On the other hand, when N, K ± , .., K n _i 
2 n , diagrammatic analysis shows that iV eff diminishes exponentially with 
the number of screens, making the alternative regime iV e fj < PjiBv) much 
more likely and resulting in low information rate BN e g (cf. Case 2, Section 4). 
In other words, a long chain of independently fluctuating media separated 
by a series of screens of sparse pinholes is detrimental to time reversal (and 
perhaps any) communication systems 
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6 Conclusions 



We have analyzed the time reversal propagation in independent-multipath 
Rayleigh-fading MIMO-channels with or without pinholes. The focus of the 
analysis is the stability condition, the multiplexing gain and the multi-screen 
pinholes effect. The main results are (i) that the stability holds when MC <C 
N e ffB for the BBFS channels and M <C N e g for the NBFF channels where 
N e g is the effective number of TRA-elements, (ii) that the optimal information 
rate Pjv under the power and noise constraints is achieved in the regime 
iVcfj A M ^> Pj (Bv) and (iii) that the effective number of TRA-elements 
is asymptotically the harmonic sum of TRA-elements and the numbers of 
pinholes on all n — 1 screens when the numbers of TRA-elements and the 
pinholes of each screen are greater than 2 n . The notion of the effective number 
of TRA elements is introduced for the first time and shown to be useful in 
analyzing stability and capacity in the presence of pinholes. 
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